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ABSTRACT 

We present the first results from a program of near-infrared spectroscopy aimed 
at studying the familiar rest-frame optical emission lines from the H II regions of 
Lyman break galaxies at 2; ~ 3. By targeting redshifts which bring the lines of interest 
into gaps between the strong 0H~ sky emission, we have been successful in detecting 
Balmer and [O III] emission lines in all five galaxies observed so far with CGS4 on 
UKIRT. The typical line fluxes are a few times 10~^^ erg s~^ cm~^, approximately one 
order of magnitude lower than the limits reached with wide-field narrow-band imaging 
surveys. 

For a Salpeter IMF and a i^o = 70 km s~^ Mpc~^, = 0.1 cosmology, the 
H/3 luminosities uncorrected for dust extinction imply star formation rates of 
20 - 270 Mq yr-i; these values are greater than those which may have been deduced 
from the ultraviolet continuum luminosities at 1500 A by factors of between ~ 0.7 
and ~ 7. Uncertainties in the shape of the reddening curve and in the intrinsic UV 
continuum slope do not allow us yet to assess accurately the level of dust extinction; 
however, on the basis of the present limited sample it appears that an extinction of 
1 — 2 magnitudes at 1500 A may be typical of Lyman break galaxies. This value is 
consistent with recent estimates of dust obscuration in star forming galaxies at z < 1 
and does not require a substantial revision of the broad picture of star formation over 
the Hubble time proposed by Madau et al. (1996). 

In four out of five cases the velocity dispersion of the emission line gas is 
o" ~ 70 km s~^, while in the fifth the line widths are nearly three times larger. Virial 
masses Myij. !=a 1 — 5 x 10^° are suggested, but both velocities and masses could 
be higher because our observations are only sensitive to the brightest cores of these 
systems where the line widths may not sample the full gravitational potential. The 
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relative redshifts of interstellar absorption, nebular emission, and Lyman a emission 
lines differ by several hundred km s"^ and suggest that large-scale outflows may be 
a common characteristic of Lyman break galaxies. The forthcoming availability of 
high resolution infrared spectrographs on large telescopes will soon allow all of these 
questions to be addressed in much greater detail. 

Subject headings: cosmology:observations — galaxies:evolution — galaxies:starburst — 
infrared:galaxies 
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1. INTRODUCTION 

Star forming galaxies at redshift z ~ 3 are now being discovered in large numbers from 
deep imaging surveys designed to reveal the break in their spectral energy distribution caused by 
the limit of the Lyman series of neutral hydrogen at 912 A (Steidel, Pettini, & Hamilton 1995, 
Steidel et al. 1996; Lowenthal et al. 1997). At z = 3 optical wavelengths sample the rest frame 
ultraviolet continuum which, being produced primarily by O and early B type stars, is a measure 
of the instantaneous star formation rate. Thus, by constructing the ultraviolet luminosity function 
of galaxies in different redshift intervals, it has become possible to sketch the global history of 
star formation (and associated metal production) in the universe over ~ 90% of the Hubble time 
(Madau, Pozzetti, & Dickinson 1998). 

In principle, the ultraviolet continuum is a more direct avenue to the SFR than other 
commonly used indicators, such as the Balmer lines — which measure the reprocessed ionizing flux 
from the most massive stars — and the far-infrared luminosity due to grain emission. Concerns 
remain, however, regarding the extent to which the observed ultraviolet luminosities are attenuated 
by the dust which is likely to be associated with the star-forming regions, given that even relatively 
small column densities of dust can modify substantially the emergent spectral energy distribution 
in the far-UV. 

The magnitude of the dust correction appropriate to Lyman break galaxies is currently the 
subject of considerable debate. On the one hand, there are claims that the ultraviolet continuum 
slope and UV-optical spectral energy distribution (SED) imply that the typical z ~ 3 galaxy 
suffers a factor of ^ 10 extinction at A = 1500 A (Meurer et al. 1997; Sawicki & Yee 1998), 
although other analyses (e.g. Trager et al. 1997) arrive at more modest estimates (factors of 2 — 6) 
from similar data. The first results from searches at sub-mm wavelengths have been interpreted 
as evidence for a significant population of dust-obscured star forming galaxies at z > 2 (Blain 
et al. 1998), but the present uncertainties in the surface density of sub-mm sources and in their 
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redshifts make such conclusions still tentative. Higher metal production rates than those implied 
by the (uncorrected) ultraviolet luminosity density at high redshift have also been advocated to 
explain the lack of evolution in the Fe abundance of rich galaxy clusters up to z ~ 0.3 (Mushotzky 
& Loewenstein 1997; Renzini 1997). 

On the other hand, dust obscuration does not appear to be a major problem for galaxies at 
z ^ 1. Observations of the infrared continuum with the Infrared Space Observatory (Flores et al. 
1998) and of the Ha emission line (Tresse & Maddox 1998; Glazebrook et al. 1998) indicate that 
the typical correction to the near-ultraviolet continuum luminosity at 2800 A is at most a factor 
of two (see also Treyer et al. 1998). Similarly, the initial results from a survey of Ha emission 
associated with damped Lyman a galaxies at z ~ 2.5 by Bechtold et al. (1998) do not point to a 
significantly higher star formation rate density than that predicted by Madau et al. (1998). 

This confusing situation is unlikely to be settled simply by improved measurements of the 
ultraviolet spectra of high redshift galaxies. As discussed below, uncertainties in the shapes of 
both the unattenuated spectral energy distribution and the reddening curve limit the accuracy 
of even the most careful estimates of ultraviolet extinction. A more profitable approach may be 
to search for Balmer emission lines, particularly Ha and H/3, associated with known galaxies at 
2; ~ 3. Since the effects of dust are less severe at rest-frame optical wavelengths, a comparison 
between the star formation rates implied by the Balmer lines and by the UV continuum should 
lead to estimates of dust obscuration which are less model dependent. An added incentive is the 
possibility of obtaining from the widths of the nebular lines direct estimates of the masses of the 
galaxies which, as argued by Steidel et al. (1996), cannot be reliably deduced from consideration 
of the ultraviolet interstellar absorption lines alone. Knowledge of the masses associated with 
Lyman break galaxies would make it possible to use their clustering properties to discriminate 
between different cosmologies (Adelberger et al. 1998; Steidel et al. 1998b). 

With these goals in mind we have begun a pilot program aimed at detecting nebular emission 
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lines from galaxies at ^; ~ 3; in this paper we report on the initial results of this work. The brighter 
objects in our spectroscopically confirmed sample of Lyman break galaxies are within reach of 
high resolution near-infrared spectrographs on 4-m telescopes and indeed we have detected Balmer 
and/or [O III] emission lines in every one of the five galaxies observed so far. With the forthcoming 
availability of similar spectrographs on 8-10 m telescopes such detections will become routine and 
may eventually resolve the vexed question of dust obscuration at high redshift. 

2. OBSERVATIONS AND DATA REDUCTION 

The near-IR sky is dominated by strong and highly variable airglow emission from the 
hydroxyl radical (e.g. Ramsay, Mountain & Geballe 1992) which normally precludes spectroscopic 
observations of faint sources. At a resolving power R ^ 2500, however, the 0H~ lines are resolved 
and a significant portion of the spectrum, particularly in the X-band, becomes accessible for work. 
With the large samples of Lyman break galaxies now available it is possible to choose objects at 
favourable redshifts such that the nebular lines of interest fall in gaps between the sky lines. The 
CGS4 spectrograph on the UK Infrared Telescope on Mauna Kea, Hawaii with a 256 x 256 InSb 
array provides both the spectral resolution and the wavelength coverage necessary for a search for 
H/3 and [O III] AA4959, 5007 at z ~ 3. 

Table 1 gives details of the galaxies observed. From our Keck LRIS spectra, which are 
reproduced in Figure 1, we measure two values of redshift for each galaxy, from Lyman a emission 
(when present) and from the interstellar absorption lines; they are listed in columns 4 and 5 
respectively (all redshifts are vacuum heliocentric). As discussed below (§5), Lyman a emission is 
generally redshifted relative to the absorption lines by up to several hundred km s~^, presumably 
reflecting large-scale outflows in the interstellar media of these galaxies (Kunth et al. 1998). 
Details of the Keck observations and of the optical and infrared imaging from which the TZ 



-7- 



magnitudes and colors in Table 1 were derived will be presented elsewhere. Here we concentrate 
on UKIRT observations in the ET-band aimed at detecting H/3 and [O III] emission lines at ^ ~ 3 
and in one case, Q0201+113 B13, Ha at z ~ 2.2.0 

The data were obtained over three observing runs in September 1996, and October and 
November 1997. In 1996 we used the 150 mm focal length camera and 150 grooves mm~^ grating 
(in second order) to give a dispersion of 200 A mm~^ (6 A pixel"^); with a 1.2 arcsec wide 
entrance slit and half-pixel stepping of the spectrum on the detector to improve the sampling, we 
achieved a resolving power R ~ 2600 corresponding to FWHM ~ 8 A in the i^-band. In the 1997 
observing runs, the improved image quality delivered by the new UKIRT tip-tilt system allowed 
us to use the longer focal length (300 mm) camera; the corresponding resolution — with the same 
grating and slit width as in 1996 — was FWHM ~ 6 A. 

At the telescope we followed standard observing techniques; a good description of the finer 
points to be taken into account when observing faint objects with UKIRT and CGS4 is given by 
Eales &: Rawlings (1993). The galaxies were acquired by blind offsets from nearby stars whose 
relative positions we had previously measured from our own CCD images with a typical accuracy 
of ±0.2 arcsec. We only used offsets stars within ~ 2 arcmin of the galaxies to make sure that 
both galaxy and star would fall on the dichroic filter used to direct the optical light to the guiding 
system. The spectra were recorded on two sets of rows on the detector separated by ~ 20 arcsec; 
the object was beam-switched between these two positions on the slit in a ABBA sequence with 
600 s integration at each position. 

The choice of 600 s as the length of an individual exposure (significantly longer than is 
normally the case when working in the infrared) was driven by the need to minimize the read-out 

®Q0201+113 B13 does not qualify as a Lyman break galaxy according to our usual photometric selection criteria. 
It is one of several galaxies in the tail of the redshift distribution of our sample, discovered in the process of exploring 
the boundaries for Lyman break galaxies in the {Un ~ G) vs. (G — TZ) color plane (Steidel et al. 1995). 
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noise of the detector relative to the sky noise in the dark regions of sky between the 0H~ lines 
where the Balmer and [O III] emission lines are expected to fall. The penalty which one pays for 
such long exposures is that some of the 0H~ lines saturate, or in any case vary sufficiently between 
successive exposures, to make it impossible to subtract them out properly. In the circumstances, 
we preferred this course of action (relative to shorter exposures), because our observing strategy 
was specifically aimed at avoiding the sky emission lines. 

After two sequences of four 600 s exposures (i.e. approximately every 4800 s) we returned 
to the offset star and realigned the optical and infrared beams (to take into account the effects 
of differential atmospheric refraction) before repeating the blind offset procedure. With total 
exposure times of between 18 000 and 40000 s (see Table 1), the final spectra are the sum of 
many such groups of four exposures. Between successive nights we stepped the grating so as to 
shift the spectra by a few pixels along the rows of the detector; this reduces the residual fixed 
pattern noise after flat-fielding. In general we summed only data obtained on photometric nights 
with ;S 1 arcsec seeing. One exception are the October 1997 observations of Q0201+113 B13 
which were obtained in patchy cloud and 1-1.5 arcsec seeing but which were nevertheless found to 
improve the flnal signal-to-noise ratio of the spectrum. 

Wavelength calibration was by reference to the emission line spectra of Ar and Kr hollow- 
cathode lamps internal to the spectrograph; we further used the sky 0H~ emission lines to 
monitor any wavelength shifts during the long series of exposures. Observations of standard 
stars, normally of spectral type AO, provided an absolute flux scale for the galaxy spectra. The 
major steps of the data reduction procedure are also well described by Eales & Rawlings (1993). 
Briefly, the UKIRT CGS4 data reduction system provides flat-fielded two-dimensional images of 
each group of four exposures where the sky background has been subtracted by subtracting the 
sum of the two frames with the object at position B from the sum of the two frames with the 
object at A. We used IRAF to wavelength calibrate and rebin these 2-D images to a linear scale; 
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we then removed the sky residuals (due mainly to temporal variations in the 0H~ lines) in each 
wavelength bin by fitting the background along the slit with a first order cubic spline. The last 
steps involved extracting the galaxy spectra at positions A and B, multiplying B by —1 (since the 
initial UKIRT pTOcessing described above produces a negative signal at position B), and adding 
the resultant two spectra. We then co-added the extractions from each group of four exposures 
to produce the final spectra. In general we found that using weighted extraction and co-adding 
algorithms made little difference to the final signal-to-noise ratio. 



2.1. Results 

Figures 2 and 3 show portions of the UKIRT spectra obtained in 1996 and 1997 respectively; in 
Figure 4 we have reproduced the inner three cross sections of the reduced two dimensional images 
from the 1996 observations. In Table 2 we list the measured redshifts and fluxes of the nebular 
emission lines covered, together with the line luminosities deduced for Hq = 70 km s~^ Mpc"-*- 
and qo = 0.1 (we adopt these values throughout the paper). From the continuum fluxes implied 
by the Kab magnitudes (see Table 1) we calculate the values of rest frame equivalent width listed 
in Table 2. It can be seen from the figures that, by targeting redshifts such that the nebular lines 
are well removed from the strongest night sky features, we have generally succeeded in detecting 
[O III] and Balmer emission lines from the H II regions of Lyman break galaxies at z ~ 3. Even so, 
it is clear that the observations are at the limit of what can be achieved with 4-m class telescopes, 
resulting in detections at levels of only a few a (see Table 2) despite the very long exposure times. 
In general [O III]A4959 is in the noise, but in each case the upper limit on its flux is consistent 
with the measured flux of A5007 which is three times stronger. In one case, Q0201-I-113 C6, we 
do detect [O III]A4959 but not A5007 because at Zem = 3.094 the latter coincides with a strong 
0H~ line at 2.04993 which is saturated in our spectra. In QOOOO— 263 D6, H/3 falls outside the 
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i^'-band window. 

Typical line fluxes are a few times 10~^^ erg s~^ cm^^, approximately one order of magnitude 
lower than the limits achieved so far in wide field surveys in the near infrared using narrow band 
imaging (e.g. Thompson, Mannucci, & Beckwith 1996). This is likely to be the reason why such 
surveys have not yet found a widespread population of star-forming galaxies at high redshift. 
Targeted searches for Ha emission at z ~ 2 in QSO fields, using narrow band filters tuned to 
the redshifts of known absorbers, have been more successful (Teplitz, Malkan, & McLean 1998; 
Mannucci et al. 1998; Bechtold et al. 1998). 



3. STAR FORMATION RATES AND DUST EXTINCTION 

In Table 3 and Figure 5 we compare the values of the star formation rate (SFR) deduced 
in the five galaxies from the measured luminosities in the H/3 line and in the far-UV continuum 
(Ao = 1500 A) respectively. At this stage we do not correct the observed fluxes for dust 
obscuration. In calculating SFRh/3 we have assumed a ratio Ha/II/3= 2.75 (Osterbrock 1989) and 
Kennicutt's (1983) calibration SFR = LHa/1.12 x 10"^^ Mq yr'^ (where Lr Q, is in erg s ) which 
is appropriate for a Salpeter IMF from M = 100 to 0.1 Mq^ For QOOOO-263 D6, where only 
[O IIIJA5007 is observed, we derive an upper limit to SFRh/3 assuming that A5007/II^ > 2. This is 
a conservative limit based on data from local star-forming regions (Terlevich et al. 1991; Stasinska 
& Leitherer 1996), as well as the three Lyman break galaxies in the present sample where both 
[O III] and H/3 are detected (see Table 2). 

Turning to the UV luminosity, it is important to realize that, although the continuum at 

Madau et al. 1998 deduced a slightly different conversion factor, SFR = Lhc/1.58 x 10"*^ Mq yr ^, using more 
recent population synthesis models. Adopting this calibration would decrease the values of SFRh/3 in column 3 of 
Table 3 by 29% . 
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1500 A is dominated by short-lived O and B stars, it does nevertheless show a modest increase with 
time in continuous star formation models. (In the case of an instantaneous burst of star formation 
the UV luminosity obviously fades with time after the starburst).^ The Bruzual & Chariot (1996) 
models predict that, for a Salpeter IMF within the above limits and solar metallicity, the UV 
luminosity at 1500 A from a region undergoing continuous star formation tends to an asymptote at 
109 years, at which point a SFR = IMq yr ^ produces i^isoo = 10^^ erg s ^ Hz ^. For comparison, 
10^ years after the onset of star formation the UV luminosity corresponding to SFR = IMq yr~^ 
is only 60% of the above value (see also Leitherer, Robert, & Heckman 1995). In the Bruzual & 
Chariot models these conversion factors do not depend sensitively on metallicity]^ 

At z ~ 3 our TZ filter (with effective wavelength A = 6850 A) samples the rest-frame 
continuum near Aq = 1710 A. Since our magnitudes are on the AB scale (AB = —48.60 — 2.5log fi^), 
the median TZ = 23.5 of the Lyman break galaxies considered here corresponds to 
fi, = 1.5 X 10^^9 erg s^^ cm~^ Hz^^; at z = 3 this in turn translates to a continuum 
luminosity Luiq = 3.1 x 10^^ erg s^^ Hz^^ . In arriving at the star formation rates listed in 
columns (4) and (5) of Table 3, we have applied small ^-corrections (generally less than 10%), 
based on the spectral slopes measured from our Keck LRIS spectra, to deduce -L1500 from the 
values of Lmo implied by the galaxy TZ magnitudes. 

There are two conclusions that can be drawn from the results in Table 3. First, in most 
cases the Balmer lines indicate star formation rates which broadly agree, to within factors of 
2 — 3, with the values obtained by applying the Bruzual & Chariot (1996) models to the observed 

^Throughout this paper we only consider 'continuous' star formation models — in which gas is converted into stars 
at a constant rate, as opposed to 'instantaneous burst' models — in which all stars are assumed to be formed at time 
t = with negligible star formation thereafter. Clearly these are two limiting cases to the real star formation history 
of Lyman break galaxies. However, in our view single burst models are less likely to apply to the galaxies observed 
here because the Balmer and [O III] emission lines would fade rapidly after such a burst and would be undetectable 
at the sensitivity of the present observations after less than 10^ years (Stasinska & Leitherer 1996). 

^An empirical determination of the metallicity dependence of the integrated UV continuum of stellar populations 
awaits the availability of the appropriate stellar libraries. 
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UV continuum luminosities. This is the level of agreement to be expected between different star 
formation indicators, even in the local universe (e.g. Meurer et al. 1995), given the systematic 
uncertainties in the calibrations and the possibility that ionized gas and early-type stars have 
different spatial distributions (e.g. Leitherer et al. 1996). 

Second, Figure 5 suggests that there may be a trend between the {G — TV) color and 
SFRh/j/SFRuv in the sense that the redder galaxies have apparently the larger values of this ratio. 
Thus Q0201+113 B13, which is among the bluest galaxies in our entire sample, has a Balmer line 
flux close to that predicted from its UV continuum luminosity, while in DSF 2237+116 C2, which 
with {G — TV) = 1.13 is close to the red limit of our selection criteria for Lyman break galaxies, 
we apparently see 4 — 7 times more H/? photons than expected. The other three galaxies are 
intermediate cases. Before reading too much into this 'trend', we must remember that we are 
dealing with only a small number of measurements from infrared spectra of low signal-to-noise 
ratio. Nevertheless, the effect is in the direction expected from dust extinction, which is more 
effective at ultraviolet wavelengths than in the optical. It is worthwhile considering, then, the 
magnitude of the UV extinction at 1500 A, ^1500 > implied by the SFRh/j/SFRuv ratios, and how 
this is related to the observed (G — TV) colors. 

3.1. Ultraviolet Extinction of Lyman Break Galaxies 

Specifically, we are going to determine the values of ^1500 which satisfy the requirement 
SFRuv = SFR^^, where SFR' are the star formation rates corrected for extinction. This would be 
a straightforward calculation were it not for the fact that we are totally ignorant of the wavelength 
dependence of dust extinction in Lyman break galaxies! Here we consider what are generally 
thought to be two limiting cases, the extinction curve determined for stars in the Small Magellanic 
Cloud (Bouchet et al. 1985; Gordon & Clayton 1998) and the 'attenuation' curve derived from 
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the integrated spectra of local star forming galaxies (Calzetti, Kinney, & Storctii-Bergmann 1994; 
Calzetti 1997a). The former rises steeply towards shorter wavelengths, while the latter is 'greyer' 
(see for example Figure 2a of Calzetti 1997a). These differences are thought to arise, at least in 
part, from different geometrical configurations of the dust and background sources of light. The 
SMC extinction curve is derived from observations of individual stars, where dust is mostly in a 
foreground screen which removes both scattering and absorption components of the extinction 
curve, whereas the integrated stellar spectrum which emerges from an extended region presumably 
includes photons scattered into the line of sight by dust mixed with the stars. Whether the 
latter situation also applies to Lyman break galaxies depends to some extent on the geometrical 
configuration of the large scale outflows which seem to be a common feature of their interstellar 
media (see §5 below). 

For the Bouchet et al. (1995) SMC reddening law with Pel's (1992) normalisation 

Aisoo = 3.44 X log (SFRjj^/SFRuv) (1) 

where A1500 is in magnitudes, whereas the Calzetti (1997a) attenuation curve gives 

^1500 = 4.87 X log (SFRh^/SFRuv). (2) 

With these conversion factors we deduce the values of ^1500 listed in Table 4. It can be readily 
seen that the uncertainties in the conversion from L1500 to SFR and in the wavelength dependence 
of the extinction combine to give a wide range of possibilities for ^1500- UV extinctions by ?a 1 — 2 
magnitudes seem typical, but we also find some unlikely solutions. The large dust corrections 
indicated for DSF 2237+116 C2 if the Calzetti reddening curve applies would imply star formation 
rates in excess of 1000 Mq yv~^ . We have retained negative (and therefore unphysical) values 
of A1500 in Table 4 because they give an indication of the inherent limitations of the analysis, 
which forces SFRyy = SFR^^. Other possibilities are that our assumption that A5007/H^ > 2 in 
QOOOO— 263 D6 is incorrect, and that in Q0201+113 B13 the Balmer lines suffer higher extinction 
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than the stellar continuum (or some of the ionizing photons escape the nebula). Note that the 
combination of the 10^ year old continuous star formation model, which gives the lower SFR for 
a given UV luminosity, and the Calzetti attenuation, which has the lower differential extinction 
between 1500 A and H/3, is the one which produces the largest values of ^1500 ■ 

3.1.1. Extinction Estimates from the Slope of the Ultraviolet Continuum 

With a larger sample of emission line measurements, it will be possible in future to obtain 
an independent estimate of the star formation rate density at high z which is less sensitive to 
dust corrections than the current rest frame UV data. For moment, however, it is of interest to 
consider whether the few Balmer line detections available can be used to 'calibrate', as it were, a 
dust extinction index based on the ultraviolet spectral energy distribution (SED), given the large 
number of UV spectra of Lyman break galaxies now being obtained (Steidel et al. 1998b). 

An index often used is /3, the slope of the UV continuum approximated by a power law of 
the form oc A^. Model spectral energy distributions (e.g. Leitherer & Heckman 1995) show 
that (3 changes relatively little with age and metallicity; for a Salpeter IMF from M = 100 to 
0.1 Mq the Bruzual & Chariot (1996) continuous star formation models have /3 increasing from 
~ —2.5 to ~ —2 over the period 10^ to 10^ years. Empirically, the measured values of /? correlate 
both with the far-infrared excess (Meurer et al. 1995) and with the reddening deduced from the 
Balmer decrement (Calzetti 1997a), although there is considerable scatter in these relations. The 
empirical starburst template spectrum derived by Calzetti (1997b) has (3 = —2.1 which would 
suggest that the youngest model SEDs may be too blue. 

At z ~ 3 our (G — TV) color provides a more reliable estimate of the continuum slope between 
1190 and 1710 A (the effective wavelengths of the two filters are 4740 and 6850 A respectively) 
than the value which could be measured directly from the Keck spectra. The reason for this is that 
the spectra are obtained through multi-object apertures which have fixed orientation on the sky; 
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differential atmospheric refraction can lead to increasing light losses with decreasing wavelength 
which may be difficult to calibrate out. The net effect is to artificially redden the spectra (by an 
unknown amount) and it is therefore dangerous to use published Keck LRIS spectra to measure (3. 
For comparison with other analyses, we have listed in column (3) of Table 5 the values of /3 implied 
by the observed (G — TV) colour for each of the five galaxies in the present sample. They were 
calculated by convolving power laws in F\ with the G and IZ filter transmission curves and folding 
in the average Lyman a forest opacity as a function of redshift according to the prescription by 
Madau (1995). 

In order to estimate the dust reddening of the UV continuum, we have determined a color 
excess E(G — TZ) = {G — TVjohs — {G — 7^)caic where the two terms on the right-hand side of the 
equality are the observed and expected values of {G — IZ). Again, we have calculated {G — 7^)caic 
by convolving the appropriately redshifted Bruzual & Chariot (1996) model SEDs with the filter 
transmission curves and allowing for the Lyman a forest opacity. The resulting values of E(G — IZ) 
for the 10'' and 10^ year model SEDs are listed in columns (4) and (7) of Table 5 respectively. 
Columns (5) and (6), and (8) and (9) of Table 5 then give the values of ^1500 implied by the 
(G — TZ) color excess, depending on whether the SMC or Calzetti extinction law is adopted. Once 
again we see that A1500 is not well constrained, the different permutations of intrinsic UV slope 
and reddening curve allowing solutions which span a 1 — 1.5 magnitude range. In this case it is 
the combination of the Calzetti curve and 10^ year model (with its very blue intrinsic slope) which 
requires the largest values of ultraviolet extinction. 

We can proceed further, however, by considering which set of model assumptions results in 
the values of ^1500 for each galaxy which are most consistent between Tables 4 and 5. Comparing 
the entries in the two tables we see that for QOOOO-263 D6 the combination of 10^ years SED 

Along any particular sight-line the Lyman a forest opacity may differ somewhat from Madau's average. However, 
the ensuing error in slope of the UV continuum is small compared with that which could result from the uncertain 
photometric calibration of the LRIS spectra. 
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(UV9) and SMC extinction law satisfies both the UV slope and the SFRh/j/SFRuv tests with 
^1500 ^ 0.5 . UV7 + SMC gives a good solution in Q0201+113 C6 for ^1500 ^ 0.7. The agreement 
is less than optimum in B2 0902+343 C6 and DSF 2237+116 C2, but even here we see that 
the UV7 SED gives answers which differ by only 0.5 magnitudes between Tables 4 and 5 
independently of the extinction curve used. Q0201+113 B13 is not useful in this context because 
its spectrum is essentially unreddened. 

It would be inappropriate to take these comparisons too literally and, for example, conclude 
that a particular extinction law is favoured or deduce an age for the star forming galaxies. The 
reasons are that our infrared spectra are of low signal-to-noise ratio and the galaxies probed are 
not sufficiently reddened to discriminate positively between the models. Furthermore, we do not 
even know at this stage that the models are internally consistent. For instance, the Bruzual & 
Chariot (1996) SED for the 10^ year old continuous star formation case may have a UV slope 
which is too blue compared to the unreddened continuum of the typical Lyman break galaxy and 
yet provide the correct calibration of the luminosity at 1500 A as a function of SFR — the slope 
and normalisation of the models are independent parameters to some extent. 

What seems to emerge from the comparison of Tables 4 and 5, however, is that the values of 
^1500 which fit best both the UV slopes and the ratios of the Balmer lines to the UV luminosity 
are in every case intermediate between the extreme possibilities allowed by the different models. 
Very high and very low values of UV extinction seem to be excluded even by the present limited 
sample. 

4. LINE WIDTHS AND MASSES OF LYMAN BREAK GALAXIES 

In column (2) of Table 6 we have listed the values of the velocity dispersion a (= FWHM/2.355) 
deduced by fitting the Balmer and [O III] emission lines with Gaussian profiles and correcting for 
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the instrumental resolution (thermal motions make a very small contribution to the line widths 
found here). In general there is good agreement between values of a measured from different 
emission lines in the same galaxy; the errors quoted in Table 6 were obtained by propagating 
the error of each Gaussian fit. We find that four out of the five galaxies have similar velocity 
dispersions, u ~ 60 — 80 km s~^, while in the fifth, DSF 2237+116 C2, the lines are apparently 
about three times wider. 

HST NICMOS images of two of the galaxies in our sample (see Table 6) show well resolved 
objects with half-light radii r ~ 0.25 arcsec (corresponding to ~ 2 kpc at z ~ 3 in the cosmology 
adopted here). For two other galaxies, we measure similar sizes from near-IR images obtained in 
0.5 arcsec seeing with NIRC on the Keck II telescope (when both are available, ground-based and 
ifS'T measurements are in good agreement). Combining these radii with the measured velocity 
dispersions we obtain virial masses M^ir w 1 - 5 X 10^° Mq (see column (5) of Table 6). Thus, 
in Lyman break galaxies at z ~ 3 we already find masses comparable to, or exceeding, that of 
the Milky Way bulge (Dwek et al. 1995). In reality, we may well be underestimating the total 
masses of these systems because at the low S/N of our observations we could easily miss broader 
components of the emission lines^ As discussed by Steidel et al. (1998a — see also Giavalisco 
et al. 1998 and Adelberger et al. 1998) the clustering properties of Lyman break galaxies, 
interpreted within the cold dark matter scenario of galaxy formation, point to dark halo masses 
Mdm ^ 10" Mq. 

If star formation takes place in rotationally supported disks, the typical a ~ 70 km s^^ implies 
rotational velocities Vrot — 120km s^^ (Rix et al. 1997). However, it has been known for a long 
time that in local starburst galaxies the nuclear emission lines do not reflect the full rotation speed 
of the galaxy (Weedman 1983). Lehnert & Heckman (1996) found that starbursts are roughly 

^^It is also conceivable that we may o«;er-estimating the virial masses if the emission line widths have significant 
contributions from gas accelerated to high velocities from supernova explosions and stellar winds. 
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confined to solid body part of the galaxy rotation curve and that the widths of the emission lines 
typically underestimate the full Vrot by a factor of 2. If a similar situation applied to galaxies at 
2; ~ 3 (a highly speculative hypothesis, admittedly) rotational velocities v^-ot ~ 200 — 250 km s~^ 
may be indicated. Note that in this picture, where the emission in the nebular lines is dominated 
by H II regions in the star forming cores of larger galaxies, it is reasonable to expect the strong, 
saturated interstellar absorption lines to be broader than H/3 and [O III], as is indeed the case 
(Steidel et al. 1996), because the absorption takes place over longer pathlengths, presumably 
through half of the galaxy. In this case both galactic rotation and the large scale 'stirring' of 
the interstellar gas by supernovae and stellar winds (sec §5 below) may be contributing to the 
strengths of the absorption lines. Heckman (1998) reached a similar conclusion for local starburst 
galaxies. 

5. REDSHIFTS OF EMISSION LINES AND LARGE SCALE MOTIONS 

Table 7 lists the relative velocities of the interstellar absorption, nebular emission, and 
Lyman a emission lines; we have assumed that the Balmer and [O III] emission lines are at 
the systemic redshift of each galaxy. The immediate conclusion is that large velocity fields are 
a common feature of Lyman break galaxies, confirming the initial hints provided by the large 
equivalent widths of the interstellar absorption lines (Steidel et al. 1996) and in agreement with 
the analyses by Lowenthal et al. (1997) and Pranx et al. (1997). In all three cases where we 
detect Lyman a emission, the line peak is redshifted by « 1000 km s~^ relative to the metal 
absorption lines. In two out of three cases (QOOOO-263 D6 and B2 0902+343 C6) H/3 and [O III] 
are at intermediate velocities; furthermore, in QOOOO— 263 D6 Lyman a emission exhibits an 
obvious P-Cygni profile, with a sharp drop on the blue side and a long tail of emission extending 
to ~ +1100 km s~^ on the relative velocity scale of Table 7 (see Figure 8 of Pettini et al. 1997). 
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These characteristics are most easily interpreted as evidence for large scale outflows with 
velocities of ^ 500 km in the interstellar media of the galaxies observed. In this picture 
Lyman a emission is suppressed by resonant scattering and the only Lyman a photons which can 
escape unabsorbed in our direction are those back-scattered from the far side of the expanding 
nebula, whereas in absorption against the stellar continuum we see the approaching part of 
the outflow. The data also show that the real situation is probably more complex than this 
simple sketch. In one case (DSF 2237+116 C2) H/? and [O III] emission are apparently at 
roughly the same velocity as the absorption lines, even though Lyman a emission is redshifted by 
~ 1000 km s^^. Most difficult to understand in the 3200 km s~^ difference between emission and 
absorption lines in Q0201+113 C6, large enough to raise the question of whether the two sets of 
lines are at all related! 

Nevertheless, all of the above features — blueshifted metal lines, redshifted Lyman a 
emission, and a variety of Lyman a profiles ranging from emission to P-Cygni and to damped 
absorption — have been observed in HST spectra of nearby H II and starburst galaxies (Kunth et al 
1998; Gonzalez Delgado et al. 1998), although lower outfiows velocities are normally involved. The 
variety of Lyman a emission-absorption profiles may refiect different stages in the interaction of 
the mechanical energy generated by the starburst with the interstellar medium of the host galaxy 
and/or different viewing angles (Giavalisco, Koratkar, &; Calzetti 1996). Irrespective of the details 
of this picture, two facts are clear. First, we see directly the process by which heavy elements can 
be distributed far from the sites of production (Heckman 1998). If outflows of ~ 500 km s~^ 
are maintained over lifetimes of ^ 1 x 10® years, the metals will travel over distances ^ 50 kpc, 
more than enough to seed the entire dark matter halo of each Lyman break galaxy. Second, just 
as the luminosity of Lyman a emission is in general poorly related to the star formation rate, its 
wavelength is not a useful measure of the galaxy systemic redshift. 
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6. DISCUSSION 

Our first detections of nebular emission lines at z ~ 3 have predictably gone only part way 
towards determining the typical extinction suffered by the ultraviolet continuum of Lyman break 
galaxies. The present sample of 5 galaxies includes one which is essentially unreddened, three 
where the continuum at 1500 A is dimmed by between ~ 0.5 and ~ 1 magnitude, and one where 
the extinction may be as much as ~ 2.5 magnitudes. The main conclusion reached in §3, however, 
is that the values of ^1500 which best fit all the available data are intermediate solutions within 
the wide range allowed by different combinations of SEDs and reddening curves. 

Recently, Dickinson (1998) has analysed the (G — TZ) colors of our entire sample of more 
than 400 spectroscopically confirmed Lyman break galaxies in the same way as discussed in §3.1.1 
above. His analysis showed that, based on the UV slope alone, the value of ^1500 to be applied 
to the population as a whole can be as little as 0.75 magnitudes or as large as 3.1 magnitudes for 
the difi^erent combinations of SEDs and extinction curves considered in §3.1.1 . The upper end of 
this range agrees with the factor of 16 proposed by Sawicki &; Yee (1998); in both cases such high 
values are derived assuming the bluest intrinsic continuum (UV7 in the notation above) and the 
greyest extinction curve (Calzetti). 

If our preliminary conclusion for the five galaxies considered here applies to the whole 
population, it would favour values of ^1500 near the middle of the range determined by Dickinson; 
thus, A1500 2i 1 — 2 may be the typical ultraviolet extinction (in magnitudes) suffered by Lyman 
break galaxies at 1500 A. If this still tentative statement holds up in the light of future infrared 

observations, several interesting consequences follow. 

First, this amount of dust would produce an extinction ^2800 — 0.5 — 1.5 magnitudes (SMC 
and Calzetti curves respectively) of the near-ultraviolet continuum of galaxies at 2; ^ 1 in the 
Canada-Prance Redshift Survey (Lilly et al. 1996). This is in good agreement with available 
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determinations of >l2800) based on the comparison between the luminosities in Ha and in the 
continuum at 2800 A (Tresse & Maddox 1998; Glazebrook 1998; see also Flores et al. 1998). 
Similarly, Buat k. Burgarella (1998) recently concluded that ^2000 — 1-2 magnitudes is typical of 
nearby starburst galaxies. Thus there is no need to assume that, as we look back to earlier epochs, 
an increasing fraction of the star formation activity takes place in highly obscured galaxies (such 
a scenario would be called for if the far-UV continuum of Lyman break galaxies were suppressed 
by more than a factor of 10; Madau et al. 1998, Guiderdoni et al. 1997). The dust corrections 
proposed here, while raising the values of the ultraviolet luminosity density by larger factors at 
z > 3 than at z < 1, still maintain the same broad picture of the cosmic star formation history 
sketched by Madau et al. (1996), with the peak in activity at an intermediate epoch between these 
two redshifts. 

(Of course the present data do not address the separate question of whether our photometric 
selection based on the Lyman break misses a significant component of the galaxy population at 
2; ~ 3, in which the the UV continuum may be too faint to be detected. However, as discussed by 
Adelberger et al. (1998), the clustering properties of the Lyman break galaxies suggest that such 
an unseen population is likely to contribute only a small fraction of the total number of luminous 
galaxies at these redshifts.) 

Second, it is of interest to ask where the metals associated with the higher, dust-corrected 
rates of star formation are to be found. Pettini et al. (1997) showed that at z ^ 2 damped 
Lyman a systems can account for the metal production expected from the uncorrected UV 
luminosity. We note, however, that the estimates of metallicity by Pettini et al. are based on 
measurements of the abundance of Zn, which traces that of the Fe-peak elements. The average 
metallicity of DLAs, Zqla^ may need to be corrected upwards by a factor of ~ 2, if O and other a 
elements are overabundant by 0.4 dex relative to Fe, as is the case in metal-poor stars of the Milky 
Way (e.g. McWilliam 1997). Therefore even when dust obscuration is taken into account, there 
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may still be a broad consistency, given the uncertainties, between these two independent measures 
of the metallicity of the universe — the metals seen in absorption towards distant QSOs and those 
produced by the star forming galaxies imaged directly. On the other hand, if future observations 
were to reveal that we have underestimated the dust correction to the UV luminosity of 2; ~ 3 
galaxies, the implication would be that DLA systems 'miss' metal-rich regions of the universe, 
either because current QSO absorption line samples are biased in favour of metal- and dust-poor 
sightlines, and/or because the metals may be retained close to their production sites. It will be 
possible to test these ideas in the near future, as wide field surveys yield new and larger samples of 
damped Lyman a systems towards fainter — and potentially more reddened — QSOs than studied 
up to now (e.g. Shaver et al. 1998). 



7. SUMMARY 

The main results of this work are as follows. 

(1) We report the first detections of Balmer and [O III] emission lines from five Lyman break 
galaxies at z ~ 3, achieved by targeting redshifts such that these transitions fall in the gaps 
between the strong 0H~ lines which dominate the near-infrared sky. The five galaxies span nearly 
the full range of our color selection criteria for Lyman break galaxies. 

(2) For a Salpeter IMF between 100 and 0.1 M© and neglecting dust extinction, the measured 
H/3 luminosities of ~ 0.8 - 10 x 10^^ erg s"^ {Hq = 70 km s'^ Mpc-^ qo = 0.1) imply star 

formation rates of between ~ 20 and ~ 270 Mq yr~^. These values are generally larger than those 
estimated from the ultraviolet continuum at 1500 A by factors of between ^ 0.7 and ~ 7. 

(3) The present sample is too small to determine reliably the value of v4i5oo, the dust 
extinction suffered by the continuum at 1500 A, given the uncertainties in the model spectral 
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energy distributions and in the wavelength dependence of the extinction curve. However, in 
general we find that the solutions which best fit both the slope and the luminosity of the UV 
continuum (relative to H/3) correspond to intermediate values of extinction within the large range 
allowed by the models. Thus, we favour A1500 2± 1 — 2 magnitudes as the values most likely to be 
representative of the whole sample of Lyman break galaxies. 

(4) The amount of dust implied would depress the near-UV continuum by only 

^2800 — 0.5 — 1.5 magnitudes, in good agreement with recent determinations of this quantity for 
CFRS galaxies at 2; < 1. The upward corrections proposed do not require a significant revision 
of the broad picture of the star formation history of the universe put together by Madau et al. 
(1996). Similarly, there is still a plausible agreement between the metal production associated with 
the star forming galaxies we see and the heavy elements detected in absorption against distant 
QSOs. 

(5) The nebular emission lines recorded are resolved in every case. It appears that the 
typical velocity dispersion is a ~ 70 km s~^, although we also find one galaxy where the line 
widths are nearly three times larger. With typical half-light radii of ~ 2 kpc, virial masses 
Myir 1 — 5 X 10-*^^ Mq are suggested. It is possible that both velocities and masses have been 
underestimated because with the limited sensitivity of the present observations we may well be 
seeing only the inner cores of the galaxies, where star formation is most in evidence. 

(6) The relative redshifts of the interstellar absorption, nebular emission, and Lyman a 
emission lines indicate that large scale outflows with velocities of at least ^ 500 km are a 
common feature of the interstellar media of Lyman break galaxies. These outflows, presumably 
driven by the mechanical energy generated in the star formation episodes, can distribute the 
products of stellar nucleosynthesis over large volumes. 

(7) Perhaps the most valuable aspect of this work is in showing that infrared spectroscopy 
of galaxies at z 2± 3 is feasible and in highlighting the wealth of information which such spectra 
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potentially offer. With the forthcoming availabiUty of medium resolution spectrographs on 
large telescopes, we expect that infrared observations will play a major role in advancing our 
understanding of the nature of Lyman break galaxies. 
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TABLE 1 
GALAXIES OBSERVED 



Name 




RA 


a 


Dec" 




^Lya 


~c 

^abs 


7^ 


(G-TZ) 


(7^ - Kab)'^ 


Exp. Time (s) 


UKIRT Obs. Run 


QOOOO-263 D6 


00 


03 


23.8 


-26 02 


49 


2.971 


2.961 


22.88 


0.45 


0.42 


21600 


Sep. 


1996 


Q0201+113 C6 


02 


03 


41.8 


+11 34 


42 




3.053 


23.90 


0.58 


0.50 


18 000 


Sep. 


1996 


Q0201+113 B13 


02 


03 


49.3 


+11 36 


11 




2.168 


23.43 


-0.06 


0.93 


39 500 


Oct. & Nov. 1997 


B2 0902+343 C6 


09 


05 


20.5 


+34 09 


08 


3.099 


3.080 


24.13 


0.45 


1.21 


21600 


Nov. 


1997 


DSF 2237+116 C2 


22 


40 


08.3 


+11 49 


05 


3.333 


3.319 


23.55 


1.13 


1.25 


33 600 


Nov. 


1997 



''J2000 coordinates 

''Redshift of Lyman a emission when present 
"^Redshift of interstellar absorption lines 

"^AU magnitudes are in the AB system; i^AB = Kg + 1.86, where Kg is the magnitude on the Vega scale 



TABLE 2 



REDSHIFTS, FLUXES, LUMINOSITIES, AND EQUIVALENT WIDTHS OF NEBULAR EMISSION LINES 



Name H/3 A4861.32 [0 III] A4958.91 [0 III] A5006.84 Ha A6562.82 







pa 
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pa 








pa 








pa 






(1) 
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(8) 
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QOOOO-263 D6 


d 


d 


d 


d 




< 2.5° 


< 2' 


< 23' 


2.966 


7.6 + 0.7 


6.1 + 0.6 


70 + 6 










Q0201+113 C6 


3.094 


2.6 + 0.6 


2.3 + 0.5 


53 + 11 


3.094 


4.6 + 0.9 


4.0 + 0.8 


95 + 20 


/ 


/ 


/ 


/ 










Q0201+113 B13 


























2.166 


6±1 


2.1 ± 0.4 


72 ± 10 


B2 0902+343 C6 


3.091 


3.0 + 1 


2.7 + 0.8 


40 + 13 


3.092 


3.3 + 1 


3.0 + 0.9 


44 + 13 


3.091 


7.7 + 1.1 


7.0 + 1 


102 + 15 










DSF 2237+116 C2 


3.316 


10 + 2 


11 + 2 


80 + 13 




< T 


< T 


< 55' 


3.318 


33 + 5 


37 + 5 


260 + 40 











































Notes. — Emission redshifts are vacuum heliocentric, la errors were measured from the error spectra shown in Figures 2 and 3. 
" Line flux in units of lO^^'^ erg s^^ cm^^. 

' Line luminosity in units of Iff*^ erg s^^ (Hq = 70 km s^^, go = O-l). 
° Rest frame equivalent width in A. 

Outside the /f-band. 
' 3f7 limit. 

^ Blended with sky OH^ emission. 



TABLE 3 

STAR FORMATION RATES UNCORRECTED FOR DUST EXTINCTION 



Name 


(G-7^) 


SFRh^ (Mq yv-y 


SFRuvT (Mq yr-i)'' 


SFRuv9 (Mq yr-i)^ 


QOOOO-263 D6 


0.45 


< 75"' 


85 


49 


Q0201+113 C6 


0.58 


55 


34 


20 


Q0201+113 B13 


-0.06 


19^ 


29 


17 


B2 0902+343 C6 


0.45 


67 


32 


19 


DSF 2237+116 C2 


1.13 


270 


65 


37 



'^Star formation rate deduced from the H/3 luminosity; typical error in the measured H/3 flux is < 30% (see Table 2). 
^Star formation rate deduced from the continuum luminosity at 1500 A adopting a 10^ year old continuous star 
formation model (see text for other assumptions). 

•^Star formation rate deduced from the continuum luminosity at 1500 A adopting a 10^ year old continuous star 

formation model. 

Assuming H/3/A5007 < 0.5. 
'^Deduced directly from Ha 



TABLE 4 

EXTINCTION AT 1500 A DEDUCED BY COMPARING BALMER LINE AND UV CONTINUUM LUMINOSITIES 



Name 


(G-7^) 


^uv7 (SMC)" 


Auvr (Calzetti)* 


^uv9 (SMC)'= 


^UV9 (Calzetti)"^ 


QOOOO-263 D6 


0.45 


< -0.19 


< -0.26 


< 0.64 


< 0.90 


Q02Q1+113 C6 


0.58 


0.72 


1.02 


1.51 


2.14 


Q0201+113 B13 


-0.06 


-0.63 


-0.89 


0.17 


0.24 


B2 0902+343 C6 


0.45 


1.10 


1.56 


1.88 


2.67 


DSF 2237+116 C2 


1.13 


2.13 


3.01 


2.97 


4.20 



''lO year old continuous star formation + SMC extinction curve. 
''lO^ year old continuous star formation + Calzetti attenuation curve. 
"^10^ year old continuous star formation + SMC extinction curve. 
''lO^ year old continuous star formation + Calzetti attenuation curve. 



TABLE 5 

EXTINCTION AT 1500 A DEDUCED FROM THE SLOPE OF THE UV CONTINUUM 



Name 


(G-n) 


(:V' 


E(G-7^)'' 


Ajjyr (SMCy 


Auyr (Calzetti)"' 


E{G-TIY 


^uv9 (SMC)/ 


^UV9 (Calzetti)^' 


QOOOO-263 D6 


0.45 


-1.1 


0.41 


0.80 


2.13 


0.26 


0.50 


1.33 


Q0201+113 C6 


0.58 


-1.2 


0.43 


0.72 


1.77 


0.22 


0.36 


0.88 


Q0201+113 B13 


-0.06 


-2.2 


0.04 


0.11 


0.23 


-0.09 


-0.27 


-0.53 


B2 0902+343 C6 


0.45 


-1.5 


0.27 


0.44 


1.08 


0.05 


0.09 


0.21 


DSF 2237+116 C2 


1.13 


-0.4 


0.71 


1.14 


2.44 


0.46 


0.75 


1.58 



"^Spectral slope of the UV continuum {Fx oc A'') after correction for Lyman a forest opacity. 

'^Excess over the color predicted by the model SED for a 10^ year old continuous star formation and Lyman a forest opacity. 
•^10^ year old continuous star formation + SMC extinction curve. 
■^10^ year old continuous star formation + Calzetti attenuation curve. 

'^Excess over the color predicted by the model SED for a 10^ year old continuous star formation and Lyman a forest opacity. 
'^lO^ year old continuous star formation + SMC extinction curve. 
^10^ year old continuous star formation + Calzetti attenuation curve. 
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Table 6. VELOCITY DISPERSIONS AND VIRIAL MASSES 



Name 


a (km s ^) 


Half-light radius (arcsec) 


Half-light radius (kpc)" 


Mvi, (IOIOMo) 


QOOOO-263 D6 


60 ± 10 


0.22*' 


1.8 


0.8 


Q0201+113 C6 


70 ±20 


0.25'' 


2.0 


1.2 


Q0201+113 B13 


85 ± 15 


0.2= 


1.5 


1.3 


B2 0902+343 C6 


55 ± 15 








DSF 2237+116 C2 


190 ± 25 


0.2= 


1.5 


5.5 



>^Ho = 70 km s-l; go = 0.1 

''Prom HSTNICMOS images (Giavalisco et al., in preparation) 
=Prom Keck II NIRC images 



Table 7. RELATIVE VELOCITIES INTERSTELLAR ABSORPTION, NEBULAR EMISSION, 

AND LYMAN a EMISSION LINES 



Name 


l^is abs (km s 1) 


Vh II (km s 1) 


^Lya (km S 1) 


QOOOO-263 D6 


-435 ± 50"^ 





+375 


Q0201+113 C6 


-3020 ± 30 







Q0201+113 B13 


+250 







B2 0902+343 06 


-800 ± 300 





+575 


DSF 2237+116 02 


+110 + 60 





+1090 



'^The errors quoted are the standard deviation of different interstellar absorption 
lines. 
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Fig. 1. — LRIS optical spectra of the five galaxies studied in this work. Each spectrum is shown 
in the rest frame at the redshift of the interstellar absorption lines; the flux scale is arbitrary. The 
positions of the most prominent interstellar lines are indicated at the top. The spectra have been 

smnntViprl witln a Vprnpl nf wiHtVi 19 A (thf aner-tra] rpsnlntinn^ fnr Hisnlnv 
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Fig. 2. — Portions of the UKIRT spectra of z ~ 3 Lyman break galaxies secured during the 
September 1996 observing run. The positions of the nebular emission lines covered are indicated. 
The short-dash line shows the la error applicable to each spectrum. The vertical long-dash lines 
mark the locations of the strongest sky OH^ emission features; although they have been subtracted 
out, large residuals can remain if the sky lines are saturated (see text). 
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Q0201 + 113 B13 z^^=2.166 
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Fig. 3. — As for Figure 2, showing the spectra recorded in the 1997 observing runs. 




Fig. 4.— Portions of the co-added two dimensional CGS4 images of QOOOO-263 D6 and Q0201+113 
C6. Each spatial increment is 1.2 arcsec along the slit; the position of each galaxy on the slit was 

adjusted so that most of the light is in the central cross-section. Most of the bright pixels are 
residuals from the subtraction of the strongest 0H~ sky emission lines, the positions of which are 
indicated in Figure 2. 
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Fig. 5. — The ratio of the star formation rates imphed by the luminosities of the H/3 hne and of the 
UV continuum respectively is plotted as a function of the observed {G — TZ) colour. Open circles: 
10^ year old continuous star formation model; filled circles: 10^ year old continuous star formation 
model. The horizontal dashed line is at SFRjj/j = SFRuv- 



